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It is shown in detail how the ground-state self-energy T,(k,uj) of the spin-unpolarized 
uniform electron gas (with the density parameter r s ) in its high-density limit r s — * 
determines: the momentum distribution n(k) through the Migdal formula, the ki- 
netic energy t from n(k), the potential energy v through the Galitskii-Migdal formula, 
the static structure factor S(q) from e = t + v by means of a Hellmann-Feynman 
functional derivative. The ring-diagram partial summation or random-phase ap- 
proximation is extensively used and the results of Macke, Gell-Mann/Brueckner, 
Daniel/Vosko, Kulik, and Kimball are summarized in a coherent manner. There 
several identities were brought to the light. 
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I. INTRODUCTION 

Although not present in the Periodic Table the homogeneous electron gas (HEG) is still an 
important and so far unsolved model system for electronic structure theory, cf. e.g. In its 
spin-unpolarized version, the HEG ground state is characterized by only one parameter r s , 
such that a sphere with the radius r s contains on average one electron Q . It determines the 
Fermi wave number as kp = l/(ar s ) in atomic units (a.u. ) with a = [V^tt)] 1 / 3 « 0.521062 
and it measures simultaneously both the interaction strength and the density such that high 
density corresponds to weak interaction and hence weak correlation [3]]. For recent papers 
this limit cf. 0, [|], If], Q, [s, 1(], H]. Usually the total ground-state energy per particle is 



on 



written as (here and in the following are wave numbers measured in units of hp and energies 
in units of k F ) 

3 i 3 olt 

e = e + e x + e c , e = ---, e x = -, e c = (ar s ) 2 [alnr s + (b + b 2x ) +0(r s )), (1.1) 

5 2 4 7r 

where e is the energy of the ideal Fermi g the exchange energy in lowest (1st) order 

(the corresponding direct term is zero, because the system is neutral), and e c is referred 
to as correlation energy. The constants a and b arise from the ring-diagram summation as 
explained in the following. Naively one should expect that in the high-density limit the 
Coulomb repulsion e 2 /r jsj can be treated as perturbation. But in the early theory of the 
HEG, Heisenberg has shown [l^|, that ordinary perturbation theory with e c = e 2 + e% + • • • 
and e n ~ (ar s ) n , where the subscript n is the perturbation order, does not apply. Namely, 
in 2nd order, there is a direct term e 2 d and an exchange term e 2x , so that e 2 = e 2 d + e 2x . 
Unfortunately the direct term e 2 d logarithmically diverges along the Fermi surface (i.e. for 
vanishing transition momenta q): e 2 A — >■ Inq for q — > [14]. This failure of perturbation 



theory has been repaired by Macke [15|] with an appropriate partial summation of higher- 
order terms e3 r ,e4 r , ■ • • (the subscript "r" means "ring diagram") up to infinite order. The 

n 

result is Eq. ([HI]) with a = (1 - ln2)/vr 2 w 0.031091 after Macke [15j and b w -0.0711 
after Gell-Mann and Brueckner [16j. The latter means, that there is a pure 2nd-order 
remainder of the ring-diagram summation in addition to the non-analyticity r 2 ln?v A 
consistent description up to this order requires to take into account all other terms of the 



same order, i.e. e 2x ~ r 2 . After Onsager, Mittag, and Stephen [l2[] it is e 2x = (ar s ) 2 b 2x with 
b 2x = (1/6) In 2 — 3£(3)/(27r) 2 m +0.02418. Thus part of the direct term b is compensated 
by the exchange term b 2x = —0.34 b r\ b + b 2x = 0.66 b. The physically plausible partial 
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summation of higher-order perturbation terms used by Macke, respectively Gell-Mann 
and Brueckner is called ring-diagram summation with its characteristic particle-hole-pair 
excitations (k — > k + q, \k\ < 1, \k + q\ > 1, q = momentum transfer), also known 
as the RPA = random phase approximation. It is the simplest approximation which 
simultaneously describes the closely related phenomena of screening effects and plasma 
oscillations as well as plasmon propagation with dispersion. (For their damping one has to 
go beyond RPA with local field corrections.) In the high-density limit, correlation "c" starts 
with RPA or ring-diagram terms "r" , symbolically written as c = r + • • • . In the following 
the screening parameter q 2 = 4ar s /n (being the interaction strength ar s times 4/7r) is used. 
In terms of this screening wave number q c the electron-gas plasma frequency (measured in 
units of kp) is uj p \ = q c j \f?> = y^Aar s / (37r). The building elements of the RPA Feynman 
diagrams are shown in Figs. la,b. The diagrams for e r and e2 X are in Figs. lc,d, middle parts. 



The non-analytical behavior of the total energy e at the high-density limit carries over to 



its kinetic and potential components, t respectively v, through the virial theorem 
to the chemical potential fi through the Seitz theorem 181 ]: 



17| and 



d 



d 1 



v = r s —e , t = -r- e , \l = - - -r s — e . 



' dr 



' dr* r 



5 1 d 



3 3 dr. 



[1.2) 



In the high-density limit, this means 
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t = t + t x + t c , t = ---, t x = 0, t c = -(ar s ) 2 [a\nr s + (a + b + b 2x )+0(r s )] , (1.3) 

3 oiir 

v = v + v K + v c , v = 0, v K — — -, v c = (ar s ) 2 [2a\nr s + (a + 26 + 26 2x ) +0(r„)] , 

4 7T 

(1.4) 
-0(r s ) 
(1.5) 



1 ar s 



7T 



a In r s + ( — -a + b + b 



"2x 



Fundamental relations between the simplest quantum-kinematical quantities (momentum 
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distribution n(k) and static structure factor S(q)) and the energy components t and v are 

OO ^ OO 

t = J d(k 3 ) n(A;)y , J d(k 3 ) n{k) = 1 , < n{k) < 1 , (1.6) 



z F = n(l~) -n(l + ), < z F = 1 - 0{r s ) , ->• 00) 



A(r s 
k 8 



+ 



00 

—Jd{q 3 )[l-S{q)} q l, S(q^0) 



2u- 



+ ■ 



pi 



[l-%-oo)]-^ + 



;i.7) 



with the static 1-body quantity n(k), the momentum distribution, and with the static 2- 
body quantity 5(g), the static structure factor (SSF). The Fourier transform of 1 — S(q) 
is 1 — g(r) with g(r) > being the pair density (PD), see Sec.V. The SSF S(q) behaves 
at transition momenta \q\ — 2 non- analytically, because there occurs a topological change 
from two overlapping to two non-overlapping Fermi spheres. This causes asymptotic Friedel 
oscillations of the PD g(r — ► 00), whereas cusp singularities S(q — > 0) ~ g odd let emerge 
non-oscillatory asymptotic terms of g(r — > 00). For the asymptotic coefficients the sum rules 

00 

A(r a ) = ~£4 0(0) and B(r a ) = 2u 2 pl g(0) with 1 - g(0) = l - J d(q 3 ) [1 - S(q)} (1.8) 



hold with the on-top PD g(0) = 1/2 - 0{r s ) [19|, |20| (for n{k)) and [19|, [2l|, [22| (for S{q)). 
g(0) is given by the normalisation of 1 — S(q) and describes short-range correlations together 
with the peculiar behavior of g(r l/^c), besides it determines the large-wave-number 
asymptotics of n(k) and S(q). 



In view of ( II. 6p and (1.7), one may ask, which peculiarities of these lowest-order quantum- 
kinematical quantities cause the non-analyticities of t and v g|. The above mentioned 
drastic changes, when switching on the Coulomb interaction, show up in the redistribution 
of the non-interacting momentum distribution rio(k) = 0(1 — k) within thin layers inside 
and outside the Fermi surface |fc| = 1 and a remaining finite discontinuity z F > (Migdal 



theorem 



24 



251] ). They show up also in the plasmon behavior of S(q) within a small 



spherical region around the origin of the reciprocal space, what causes an inflexion point 
?infl, S'infl ~ ^ P i- All these reconstructions describe the long-range correlation (screen- 
ing, collective mode called plasmon), characteristic for the Coulomb interaction. They 
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are treated in lowest order again by ring-diagram summations with the replacements 
n 2d{k) — Ti r (k) and Sid(q) — > S T (q). (Note, that t 2 d arises from n 2 d(k), but t> 2( j from 
Sid(q)-) How these replacements lead to the r 2 s \ar s terms of t and v is shown in 
However, because the redistributions take place essentially only in the mentioned sensitive 
regions \\k\ — 1 1 <C g c and \q\ <C g c the r-terms approach the original d-terms far off these 
regions. Thus, to be consistent up to this order, the corresponding x-terms must be taken 
into account: n c (k) = n r {k) + n 2x (k) + ■ • • , S c (q) = S T (q) + S\ x (q) + ■ • • . These x-terms 



compensate part of the r-terms, similarly as this is the case for e c , t c 



/i c with the 



ratios of b 2x to b, a + b, a/2 + b, —a/3 + b being -0.34, -0.6, -0.43, -0.3, respectively. The 
ring-diagram summation for n(k) has been developed in 27|, |28(], an analytical extrapola- 



tion is given in [30J, for t 
S(q) has been done in 132 



re spin-polarized case see 3jJ. The ring-diagram summation for 



331 ] . In 34| both n(k) and S(q) are considered on the same footing. 



Another relevant ground state property is Dyson's self-energy H(k,u), a dynamical 1-body 
quantity. Its on-shell value (as a function of r s ) is related to the chemical-potential shift 
through the Hugenholtz-van Hove (Luttinger-Ward) theorem \i — /x = S(l,/i) [35(. Besides 
this, with the off-shell self-energy u) also the 1-body Green's function 

1 



G(k,u) 



u-t{k) -E(k,u) 



with H(k,u}) — > sign(A; — 1) i5 for r s 







;i.9) 



is known, from which follow the momentum distribution (Migdal formula 24J ) 



oo oo 

n{k) = J ^ e iuS G(k,u;) rx J d(k 3 ) n(k) = 1, t = J d(k 3 ) n(k) y 

C+ 



;i.io) 



and the potential energy (Galitskii-Migdal formula 37|], for its use in total-energy calcula- 



tions cf. 



38] and refs. therein) 



r = - / d{k 3 



du 
2tx{ 



e iu,5 G{k,cu)E(k,u) 



Xll) 



with 5^0 and C + means the closing of the contour in the upper complex tu-plane. 



Note that E(/c,o;) and thus also G(k,uj), n(k) and t, as well as v are functionals of t(k) 
and v (q). Supposed the ground state energy e is available as such a functional, then the 



(generalized) Hellmann-Feynman theorems [1 



-<*> = Tim ■ S(q) ~ 1 = m m (L12) 



hold, cf. App. A. These are equivalent writings of expressions given in 27|, |34| for n(k) 



and in 
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32] for S(q). The quantities e, n(k), t, S(q), and v result as functions of r s 
from the replacements t(k) — ► fc 2 /2 and i>(q) — > q 2 /q 2 - These fundamental relations permit 
the following procedure (see Fig. 2): If H(k,u) is available as a functional of £(fe) and 
v(q) from perturbation theory or otherwise, then n(k) can be calculated with the Migdal 
formula (11.101) . Therefrom follows t with (11.61) and v with the Galitskii-Migdal formula 
(II. lip . Finally from their sum t + v — e and the functional derivative (j!.12p the SSF S(q) 
results (and n(k) may be checked once more for consistency). So the dynamical 1-body 
quantity T,(k,u>) provides the static 2-body quantity S(q). 

Whereas in Q| the on-shell self-energy £(l,/i) has been studied, here the off-shell self-energy 
E(/c, u) is considered. The problem in 9| was, to find out the correct diagrammatic sum 
for £(l,/i) on the rhs of the Luttinger-Ward theorem, which makes it an identity in the 
high-density limit r s — > 0. The answer: the GW approximation with W= RPA and G = 
an appropriately renormalized particle-hole line yield the correct r s -behavior, which agrees 
with the lhs as it follows from ^i{r s — > 0), cf. (11.51) . Here it is shown in detail how the 
off-shell self-energy H(k,uj) for the model case r s — > yields n(k), t, v, and S(q) step by 
step according to the procedure of Fig. 2. So it is shown how n(k) and S(q), which have 
their common origin in the 2-body density matrix, are indirectly linked mutually through 
the self-energy &). 

Perturbation theory for H(k,uj) means E = S x + S c with S c = S2 + S3 + ■ • • , where 
S 2 = S 2 d + S 2x - The divergence of E 2 d is corrected by the ring-diagram summation (RPA): 
£ 2 d + ' ' ' ~~ > ^rj s° S c = E r + E 2x + ■ • • . The building elements of the RPA Feynman 
diagrams (cf. Figs, la and lb) are the Coulomb repulsion 7r 2 v(q) — > n 2 q 2 /q 2 with the 
coupling constant g 2 = 4ar s /7r and the one-body Green's function of free electrons with 

'10. (1.13) 
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From Go(k,u) follows the particle-hole propagator Q(q,r)) in RPA according to 

d 3 k 



Q{q,v) = -J ^ J ^G Q {k,u)G Q {\k + q\,u + ri) 



(1.14) 



with the result 

Q{q,v) = 



d 3 k 
Aix 



+ 



_t(k + q) - t(k) -rj — id t(k + q) - t{k) + 77 - iS 



;i.i5) 



|fe|<i, \k+q\>\ 

The denominators contain the excitation energy to create a hole with k inside the Fermi 
sphere and a particle with k + q outside the Fermi sphere. Thus, Q(q,rj) is a functional 
of t(k) with the functional derivative flA.5j) and with R(q, u) = Q(q, iqu) defining a real 
function, see ( IB.lj) . In the following the self-energy contributions S x (Sec. II), S r (Sec. Ill), 
S 2x (Sec. VI) are explicitly given as they result from the diagram rules and it is derived, what 
follows from them: n(k), t, v, e, S(q) . For (here not considered) issues as S(q,ui), e(k,uj), 



GW-approximation etc. cf. e.g. 



38 



39 



40 



41 



42] and refs. therein. Sec. V deals with 



the short-range correlation following from S(q), cf. ( 11.81) . Although the high-density spin- 
unpolarized HEG is only a marginal corner in the complex field of electron correlation, its 
ring-diagram summation gives deeper insight through rigorous theorems, how energies and 
low-order quantum-kinematics are functionally related, what should be of a more general 
interest. 



II. FIRST ORDER 

The lst-order direct term vanishes because of the neutralizing positive background. So the 
expansion of E(A;, u) starts with the lst-order exchange term 

E x (k 7 [v(q)]) = - J ^^ V (q)=-1. J d 3 q v(q). (2.1) 
\k+q\<i \k+q\<i 

Its peculiarity is, that it does not depend on uj, what makes n x (k) vanishing identically, 
because of Eq. fll.lOp with G ~ Gq^Go. This is in agreement with t x = as a consequence 
of the virial theorem (|1.2p . One may therefore conclude, that all the features of n(k) start 
with the 2nd order. But this is not true. The peculiar redistribution of the non-interacting 
momentum distribution no (A;) = 6(1 — k) due to the Coulomb repulsion makes the diconti- 
nuity jump zp = n(l~) — n(l + ) to deviate from its non-interacting value of 1 in 1st order: 
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z-p = 1 — 0.18 r s + • • • . But the corresponding kinetic energy starts with t c ~ r 2 lnr s . - The 
potential energy t> x follows again as a functional of v(q) from the Galitskii-Migdal formula 
(EI]): 



^xK<?)] 



3 

8^ 



d 3 k S x (fc, [v(q)}) 



■57T 



y ci 3 fc y d 3 q v(q) . (2.2) 



|fc|<l |fc]<i,|fc+qf|<i 
This is, because of t x = 0, also the total energy in this order: e x [v (q)] = v x [v(q) 
functional derivative (11.121) yields 



Thus the 





"3g 




(1)1 


y d 3 A; = 


2 2 


2 





e(2-g) + 6(g-2). (2.3) 



|fc|<i,|fc+q|<i 

The integral arises (for q < 2) from two overlapping Fermi spheres. For q > 2 they do not 
overlap, thus a drastical change of the topology occurs, when passing q = 2. After Fourier 
transformation, see (15. ip . this non-analyticity causes the asymptotic Friedel oscillations of 
the non-interacting PD go(r — > oo) — 1 ~ cos2r, sin2r. Another non-analyticity, namely 
the cusp singularities S (q — » 0) ~ q, q 3 make the non-oscillatory terms of g (r — > oo). The 
above integral J d 3 k is just the volume of two calottes with the height h = 1 — q/2, hence 
J d 3 k = 2 • ^h 2 (3 — h). If in Eqs. (12.ip and (12.2p the general interaction line v(q) is replaced 
by the Coulomb interaction q 2 /q 2 , then 



k. 



q 2 



9£ 
4 



1 + 



k 2 



2k 



In 



fe + 1 



jfe 



and f x 



16 



9 C 



(2.4) 



turn out. Besides E x (l) = -q 2 c /A = (4/3)e x . Note e x = f x in agreement with the virial 
theorem (11.21) . This shows how the procedure of Fig. 2 works in lowest order. 



III. SECOND ORDER: THE DIRECT TERM S 2d AND ITS RPA CORRECTION 

A. The self-energy H T (k,u>) 

In 2nd order there is a direct term (d) and an exchange term (x), see Figs. lc,d, left: 
^(k, u) = S2d(fc,w) + E2x(fc,w). The direct term diverges along the Fermi surface, i.e. for 
vanishing transition momenta go — > 0. This flaw is repaired by the ring-diagram summation 
(Fig. lc, left) with the result 

<?><?o 
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Note, that this defines a functional of v(q) and note also 

V M > i for 

l + v(q)Q(q,rj) q 2 + q 2 Q(q,v) Q 2 ' 

where the screening (or "Yukawa") term q 2 Q(q,f]) in the denominator makes the bare 

Coulomb repulsion renormalized and removes the above mentioned divergence of S 2 d(fc, w)- 

Go(k, to) and Q(q,rj) are functionals of t(k), see (11. 13ft . ( 11. 15ft . Use of Eq. (11.131) leads to 

dr] v 2 (q)Q(q,r)) 



v(q)Q(q,v) 



&(\k + q\ 



+ 



G(l-\k + q\) 



uj + 77 - \k 2 - q ■ (k + \q) + \5 LU + ri-lk 2 -q-(k + lq)- 15 _ 
In the following it is shown, how n r (k) and t T result from this E r (A;,u;). 



.(3.1) 



B. How n T (k) results from Y< r (k,ui) and t T from n r (k) 



According to the Migdal formula (ll.lOp . n(k) follows from G(k,u). In RPA it is 

G(k, uj) = G {k, uj) + G (k, uj)E r (k, uj)G (k, u) + ■ • • (3.2) 



The first term yields the momentum distribution of the ideal Fermi gas 

duj 



n (k) 



2vri 



e(jfe-i) e(i -k) 



u-hk 2 + i5 u-hk 2 -i6 



0(1 - k) 



(3.3) 



being for small r s additively corrected by the RPA expression n T (k), which follows from the 
second term of Eq. (13.21) according to 



/duj 
—G (k,uj)J: r (k } uj)Go(k,uj) . 
2m 



Using £ r (&, uj) of Eq. (13. ip it is 



1 



drj v 2 (q)Q(q,r]) 



nr(k) = — I d 3 q I Pr /, V 7 w f( k ^^)- 



(3.4) 



(3.5) 



57T j j 2vri l + v(q)Q(q,T]) 
The case "2d" appears, when the "Yukawa" term v(q)Q(q,rj) in the denominator is deleted 
('descreening'). The contour integrations are comprised [using Q(k — 1)0(1 — k) — 0] in 



f(k,q,v) 



duj 

x 



e(k-i) 



+ 



0(1 - k) 



(uj - \k 2 + \b x )(u - \k 2 + i5 2 ) (uj - \k 2 - \b x )(uj - \k 2 - i5 2 ) 
e(\k + q\-l) e{i-\k + q \) 



+ 



ijj 



+ T] - \k 2 - q(k + \q) +16 uj + r] - \k 2 - q(k + \q) - i5 



(3.6) 
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Only the combinations 0(fc — 1)0(1 — \k + q\) with a pole in the upper u-p\ane and 0(1 — 
k)Q(\k + q\ — 1) with a pole in the lower plane contribute. Next the contour integration 
along the real cu-axis is closed by half-circles in the upper, respectively lower plane. Note 

that J diojiu) — u±) = ±27ri, where Im uj± ^ 0. The result is 

c± 

f(k, q, V) = T 77— i x 12 3.7 

with Q(k,q) = ±1 for |fe| ^ 1, \k + q\ ^ 1 and otherwise. This provides with the 
particle-number conservation, because of 

f f 1 d 3 k f 

J - J L- g ( fc + H r ° « /^-W = 0- (3-8) 
\fe|>i,|fe+qf|<i |fc|<i,|fc+qr|>i / 

Here, with the replacement A; — > — k — q, the denominator of the second integral transforms 
to [r] + q{k + \q)} 2 - The second integral is identical to the first one. This is because of the 
property Q(q, —rj) = Q(q, rj) in the term in front of f(k, q, r/) in Eq. (13. 5ft . 

Next it is shown, how t r results from n r (k), starting with Eq. ( 13.5H . With the identity (A. 5) 
it can be written as 



n r (k) 



~h I d3q I I^X> 1)n+v+1(<?) Qn ^ 

n=l 



St(k) 



= __L f^f^ftir^^)^. (3.9) 

An J H J 2ni ^ n + 1 yHI 5t(k) v ; 

Next t r = 3/(47r) j d 3 k n r (k) t(k) is combined with the identity (1A.6I) . It results in 

t ' = rJ d!,q S ^ ( - l)a+1 ^ a+1{q)QM ^ ■ (310) 

' ! 71=1 

This is the contribution of the ring-diagram summation to the kinetic energy. For the total 
energy contribution e r the potential energy v r is needed. 
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C. How v T results from E r (fc, uS) and S T (q) from e r 

Equations (11.91) and (13.11) inserted into the Galitskii-Migdal formula (11.111) and the 
oj— integration performed yields 



16n J J 2ni 1 + v (q)Q(q, rj) 



3 d 3 k 



0(|fe + q\ — 1)0(1 - fc) 0(1- |fe + g|)e(jfc- 1 



+ 



77 — q(fe + iqr) -7] + q(k + \q) + \5 



(3-11) 



With the definition (I1.15P of Q(q,rj) it results 



167r 7 7 27ri 1 + v(q)Q(q, 77) 
The power-series expansion (its 1st term is i> 2 d) 

3 



/A I ^ L Y.{-l) n+1 v n+ \q)Q n+ \q,r 1 ) 



/ rf3g / ^Ki + ^)Q(?^))-^)Q(g^)] • (3.13) 



makes it better comparable with Eq. (I3.10p for £ r . Their sum yields [with — ^pj + 1 = 
the well-known RPA expression for the total energy (after Macke) 

3 

e r = 

167T 

(The 1st term of the power-series expansion gives e2d-) Note that r s de T /dr s agrees with v r 
of Eq. ( 13. 121) : virial theorem (II. 21) . By means of functional derivatives [see Appendix A and 
Eqs. ( 11.121) ] follow the SSF S T (q) and the momentum distribution n r (k). 

Indeed, S r (q) results from e v [t(k), v(q)} as 

<fo(qr) 7 2tti l + v(qf)Q(g,f7) 

cf. Fig. lc, right. If this is multiplied by v(q)/(16ir) and integrated f d 3 q [according to 
(1.7)], then v x turns out, as it should. 



The analog procedure for n T (k) is 

„ (1A - Air 5e * 1 ( rfin d7] ^(^)^(g^) S Q(q>V) 6 Q(^V) J_f(U nrl \ 

l[) 3 5t(k)~ 4 J Q 2ni l + v{q)Q(q,T)) 5t(k) ' 5t(k) " 2n J[,q,V) 



(3.15) 
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in agreement with Eq. (13.51) . as it should. 



In the 

Macke 



bl 



15fl 



owing it is shown, how Eqs. ( 13. 12113. 151) rea 
Gellmann/Brueckner flo] ]. Daniel/Vosko [27] 



ly yield the high-density results of 
Kulik y, and Kimball [32] . 



D. How the complex Q{q,rj) is replaced by the real R(q,u) and 

how v r , S r (q), and n T (k) behave for r s — > 



The replacement 77 — > igw turns the contour integration in the RPA expressions (13.121 
13 . 1 5[) into one along the real axis. Besides it is a useful trick to introduce the velocity u 
instead of the frequency rj. With the replacement v(q) —>■ q^/q 2 they take the form (note 
q d 3 q = 27rq 2 d(q 2 )): 



00 00 



(2) 



^jduj d{q 2 ) q 






00 00 

it J J q 2 + q 2 R(q,u) 





(3.16) 
(3.17) 
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oo 

^2 



(3) SM = -^qfdu K[ q ^ , (3.18) 

7T J q 2 + q 2 R(q J u) 
o 

oo oo +1 

(4) n,(k) = £fdu[ Jfa* Ri i'l I di „. e <^;« „, (3.19) 



in 



' q 2 + q 2 R{q, u) J ^ q 2 [iu - {k( + \q)} 2 



-1 

with Q(k,q,e k e q ) = Q(k,q), see (13.71) . In the following, these four RPA quantities are 
discussed in detail in the high-density limit r s — > 0. 

(1) Let us first consider the total energy e r of Eq. (3.16). The power expansion leads in 
lowest order to 

oo oo oo 
ql qo 

For the momentum transfer function I(q), to be referred to as Macke function, see ( 1C.11I) . 
Its property I(q — > 0) ~ q makes to diverge for qo^O. Vice versa this flaw of the 2nd- 
order perturbation theory is rectified by the RPA summation (3.16). How this e r behaves 
'or small r s with the result (jl.ip has been shown by Macke [15] and Gell-Mann/Brueckner 

3- 

(2) In the following their method is applied to the potential energy v T of Eq. (3.17). 
Deleting the term q 2 R(q, u) in the denominator cancels the RPA partial summation and the 
diverging expression 



oo oo 



3 4 [ . [ m 2, R 2 M 2-3qi fdq 
V, d = --q e J duj d(q )—j r - = - — J -I(q) (3.21) 

ql qo 

results. It remains to show, how to extract from Eq. (3.17) the constants C\ and c<i of 
v T = (ar s ) 2 [ci lnr s + c 2 + 0(r s )} to be compared with v c = v r + O(r^) of Eq. (11. 4p . Whereas 
Ci follows from V2d, ci results from the peculiar behavior of v T for small transition momenta 
q. Therefore one can approximate R(q,u) ~ Rq(u) + ■ ■ ■ and restrict q to q < qi, where qi 



14 



is a small (non-vanishing) momentum: 



oo 9i 




|^ c 4 y du Rl(u){Hql + tiM")] - WcM^}} 
o 

oo 

qtj du R 2 (u){[lnql + O(r s )]-H q ^(u)}} . 



3 

8^ 



(3.22) 



With the constants a, b' T defined in Appendix B it is 



4a 



v » = (ar s ) 2 [2alnr s + 2(aln h b[ - alng 2 )] + 0{r 

71 



si i 



(3.23) 



thus c\ = 2a. To find also c 2 the difference Av 2d = v 2d — v® d between the correct 2nd-order 
term of Eq. (13.221) and the first term in the expansion of v®, namely 



<n 



-2d 



[oir s 



2tt 4 J q 



dq 



TV 



-=r2 ■ 4i\q—a 



(3.24) 



has to be considered (exploiting the trick of Gell-Mann/Brueckner mutatis mutandi): 

91 



Av 2d 



or. 



2tt 4 



.90 



9() 



ar. s 



dq 

2vr 4 y ^2 



■57T 4 a 



1(g) - —-0(1 - g ) 

3 9 



91 1 

+ \+0(r 3 s ) . (3.25) 



90 



The first integral does no longer diverge for q — > 0, therefore one can set go = 0. Besides 



Av 2d = {ar s ) 2 2{b 2d + a\nq 2 1 ) + 0{r 3 s ) 



(3.26) 



shows [for &2d see (B.8)], that for r s — > the sum f° + Af 2( j does not depend on the arbitrary 
cut-off q\. 



v r = {ar s ) 2 2 



4a 

a In r s + (a In h 6 r + 6 2d ) 

7T 



0(r s 3 ) . 



(3.27) 



This has to be compared with Eq. (11.41) . Indeed the constant direct term beyond lnr s yields 



4a 

a + 2b = 2(aln h b' T + b 2d ) 

7T 



(3.28) 
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which defines b. Its value b « —0.0711 agrees with what is given in [l(|, as it should. - The 
appearence of a 2nd-order term ~ r 2 means that - to be consistent - all other terms of the 
same order contribute to the small-r s behavior of e. There is only one such term, namely 
V2 X , as treated in Sec. IV. 



(3) Next the static structure factor of Eq. (3.18) is considered. In the lowest order 
(r— >ld) it is - again with the Macke function I(q) of Eq. (IC.11I) and with u p \ = q c /V3 



uj 



2 



uj 2 uj 2 
S ld (q - 0) = -3(1- In 2 )^ + 0(l/q 3 ) , S ld (q -+ oo) = -2-f + 0(1/ q 6 ) . 

At q = 2, the non-interacting value So (2) = 1 and its discontinuity jump AS '(2) = 3/4, 
which arises from the Fermi edge, are reduced by [with 1(2), AJ"(2) from Q, Eq. (C.2)] 

5 ld (2) = -(13-16lB2)i«J, ASf d (2) = - ^ oft . (3.30) 

Note, that Si x (g) compensates part (ca. half) of Sid(g), cf. Eq. (14.71) . and note that S(q) = 
S (q) + Sid(q) + Si x (g) + • • • decreases with increasing r s for a given value of g. Whereas for 
g> g c the perturbative treatment S r (g) = Sid(g) + 0(r 2 ) holds (screening effects are not so 
important for large momentum transfers q) , for q <C g c there is a big difference between the 
'bare' Sid(g), which behaves unphysically because of I(q — > 0) ~ g, cf. Eq. (1C.11I) . and its 
renormalized counterpart S r (V), where the ring-diagram summation ameliorates the above 
mentioned flaw of Sid(g), cf. [6j, Fig. 3. For small g, the approximation R(q, u) = R (u) + - ■ ■ 
is sufficient. So S r (g) ~ —q c L(q/q c ) with 

oo 

= -y du , ° l n , r> L(y -> 0) = -y - — + — -y 4 + • • ■ , (3.31) 



7T 



y 2 + i2o(u) 4 y 2 * 20 



to be referred to as Kimball function, for its properties cf. [6j. As a consequence, the 
expression 



*(»<«0 = -4« + ^-2oa--' (3 ' 32) 



3 , q 2 _ 3_(f 
; pl 20^, 

(i) eliminates the divergence of S% d (q — > 0) and (ii) simultaneously replaces in the sum 
So(g) + S r (g) the linear term of So(g) with a quadratic one, which is in agreement with the 



plasmon sum rule 



43|, 



44j. Higher-order terms arising from the difference R(q, u) — Rq{u) 
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q 2 Ri(u) + • • • and from local field corrections beyond RPA have to kill also the cubic term 
— g 3 /16 of So(q), to substitute the coefficient of the term ~ gV^pi of S T (q) correspondingly, 



46| 



and to add a term ~ q 5 , such that 145 

2 4 5 

q q q 



V ^pi "•'pi 

A direct consequence of these replacements is the appearance of an inflexion point, which 



n 1 I I 3 
Sa < U, — < s 4 < 

' 12 1 1 20 



(3.33) 



moves for r s — ► towards the origin with g infl , ~ lu p \ = ^Aar s / (3ir) . S r (q) of (13. lSj) 
realizes the smooth transition from the l/r s )-scaling "far off" the origin [reasonably 

approximated by I(q) of (13.291) ] to the k-pq c (~ 1/v^J -scaling near the origin [reasonably 
approximated by the Kimball-function L(q/q c ) of (13.311) ]. This transition causes the 
non-analyticity v r ~ r 2 In r s 6J]. 



(4) Finally the momentum distribution (3.19) is considered. The expression 
q 2 f(k, q, iqu), see (13. 7p . is developed in the following way 



1 











u 



[i« - (*C + W 



u 



+ KK + \)q? du u + i(k( + \q) du u 2 + (k( + \q) 2 



UJ 



This allows to write Eq. (3.19) as 

n r (k) 

and 



(47T/3) 



pl F T (k) with F T (k) = I du 



d 3 q R(q, 



u 



q 3 q 2 + q 2 R(q, u) 



(3.34) 



F(k,q,u) (3.35) 



+i 



rn 9 I At e(fc,g,Q 



(3.36) 



Next the (^-integration has to be performed. The boundary conditions k ^ 1 and |fc + </| ^ 1 
mean k 2 + q 2 + 2kq( ^ 1 or C ;=s Co with Co = (1 — k 2 — q 2 )/ (2kq). So we have 



F(k > 1, q, u) = - 



i H + k 



k{K + k) 2 + u 2 



Co 



, F(fc<l,g,u) 



1 K + k 



k (k( + \q) 2 + u 2 



+i 



Co 



. (3.37) 



-1 V""s ■ 2 K 

For k > 1 the function -F(fc, g, w) is non-zero only within a certain stripe of the k — q— plane 
(see Fig. 3), namely k — 1 < q < k+1 with £ = — 1 • • • Co (area I). For k < 1 the corresponding 
areas are the triangle 1 — k < q < 1 + k with C = Co • • ' + 1 (area II) and the stripe q > 1 + k 
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with £ = — 1 • • • + 1 (area III). With k( + \q = ^p- (in area I and II) it follows 



Fi(k > l,q,u) = +~ 



F n (fc < l,g,w) 



2q 



k-\q 



|? + k 



2q 
l-k 2 

2g 



F m (/c < l,g,u) 



A- 



• 1— fe 2 \2 
■ 2q 1 

\q-k 



+ u 2 {\q + k) 2 + u 2 



\ q -kf + u 2 a q + ky + u 2 



for {k, q) in I , 

for (k, q) in II , 

for (k, q) in III . (3.38) 



This together with (13.351) is the momentum distribution in the ring-diagram summation 



27 



281 ] . see also 



301 ] . For k = only the last line for III contributes with Fm(0,q,u) 



—S(q 2 — u 2 )/(q 2 + u 2 ) 2 giving -F r (0) = —4.112335 + 1.35595 r s + • ■ ■ . If in the denominator 
of Eq. (I3.35P the "Yukawa" -term q 2 R(q,u) is deleted, then 

,4 



oo oo 



n 2d (k) 



pl 



(4tt/3) : 



;F 2d (k) , F 2d (k) = 4tt 



7 



ciii i?(g, it) F(k, q, u) 



(3.39) 



arises with F 2d {k > 1) = Fj(k) and F 2d {k < 1) = Fu(k) + Fui(k). [Note that the definition 



of F 2d (k) with F 2d (k ^ 1) ^ differs from what is used in [4|, |30j, where F(k) > for all 
k ^ 1]. The function (with u ~ temp the u-integration is replaced by an angular integration 

3) 



k+1 vr/2 

F 2d (k > 1) = +^ m [dp [R(q 



1 



F 2d (k < 1) 



k-l 

l+k vr/2 

4tt r dq 
k J q 3 

l-k 

oo tt/2 



^— tan y>) - R(q, (k - -q) tan tp)] 



l-k 2 1 
[R(q, ^ tan v?) - R(q, (k+-q) tan <£>)] 



~ J ^ J d<p[R{q,{~q-k)t&n<p) - R{q,(±q + k)tim<p)] (3.40) 



i+fc 



possesses the properties (see also [J, |29[ and note that F 2x (k) is of the same order, cf. Eq. 



F 2d (k 0) = - (4.112335 + 8.984 A: 2 H ) , F 2d (A; -> oo) = X - ^ ^ + 

F M (^l-) = ± ^(l-.n2)^ + .... 



(3.41) 



Consequently, n 2d (k) approximates n T (k) for k <C 1 and ^> 1 (far off the Fermi surface, 
where the screening effect is not so important), but it diverges near the Fermi surface 
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square-inversely as ±1/(A; — l) 2 for k ^ 1. Vice versa, this divergence is removed through 
the ring-diagram partial summation with the replacement q 3 —>■ q[q 2 + q 2 R(q,u)} in Eq. 
(I3.39p . For k near the Fermi edge, it holds [6|] 



F r (k 



2tt 



G 



life -II 



for l$k$l±y/q c 



q 2 k 2 \ q, 

with G(x) being the Kulik function (IB. 151) . Thus the discontinuity jump for r s — > 
described by (the higher-order terms are different for outside/inside the Fermi surface) 



(3.42) 
■> is 



nJl z 



10. 



pi 



2ii 



G(0) + 



10 



pi 



2tt 



G(0) + 



±0.088519 r s + • • • . (3.43) 



(4vr/3) 2 q 2 " ' (4vr/3) 2 3 

Comparison of (13.411) with (13.421) shows that the divergence at A: — ^ 1 ± is eleminated and 
replaced by the non-analytical behavior of G(x —>■ 0), see ( IB. 161) . Near the Fermi surface the 
distribution is "symmetrical" and shows a logarithmical "snuggling" with infinite slopes: 



n v (k 



nr (l±) ± ^Ei [ y^L + 3 ) \ k _ i| i n \ k _ i| + Q(k - 1) 



F T (k) of (13.351) realizes the smooth transition from the fcp(~ l/r s )-scaling "far off" the 
Fermi surface [reasonably approximated by F 2 d(k) of ( 13. 39ft ] to the /cf9c(~ l/\/r s )-scaling 
near the Fermi surface [reasonably approximated by the Kulik-function G(\k — l\/q c ), see 



(I3.42p ], This transition causes the non-analyticity of t T ~ r 2 lnr s , |6]. - From (|3. l3§ follows 
z-p = 1 — 0.177038 r s + ■ • • , what is in agreement with the Luttinger formula [25(, which 
relates the quasi-particle weight z-p directly to the self-energy Ti(k,uj): 



1 



z F 



£'(l,/x) = Re 



1 - S'(l,/i) ' ^'^ y dio 
Indeed, with the ring-diagram approximation (13. ip it becomes 



(3.44) 



s:(i, 



i 

8^ 



d 3 q 



drj v 2 (q)Q(q,rj) d 



1 



27ri 1 + v(q)Q(q, rj) drj rj - q ■ (e + \q) ± 15 



\e + q\ ^ 1 . 

(3.45) 



For r. 



it behaves as 
1 



47] 



ski 



7T 2 



1 



arctan — h 
Ro(u) u 



-0.177038 r s + 



(3.46) 



in agreement with (I3.43P and (IB. 17[) . (In 48j, z-p = 1 — 0.12 r s + - ■ ■ is claimed, instead of the 
RPA figure 0.18.) The linear behavior of zp(r s ) is an example for how higher-order partial 
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summation may create lower-order terms. Calculations beyond RPA with (13.441) have been 



done in 



49J. Calculations of z-p for r s < 55 have been done in [39j.- The strength of the 



correlation tail, i.e. the relative number of particles [with k > 1 and using (IB. 191) ] is [28] 



N r = J d(k 3 ) n r {k) = u. 
l 



3 (3/2) 
pl 2tt 2 



5/2 



du 



u In • 



u 



+ ■ 



0.05383 r 3/2 + 



(3.47) 

281 ] . The M-integral is 1.06252 . How does n 2x (k) change the above results for zp and iV ? 



IV. SECOND ORDER: THE EXCHANGE TERM £ 2x AND ITS 

CONSEQUENCES 



As already mentioned above and stressed by Geldart [49j, for a consistent small-r s descrip- 
tion up to terms ~ r^lnr s and ~ r 2 the exchange terms t> 2x , e 2x , n 2x (k), Si x (q) are needed. 
Here it is shown, how they arise from E 2x (/c,u;), cf. Fig. Id, left. 



The self-energy in the second order of exchange is 
q\ f d 3 q l d 3 q 2 f di] x dri 2 



E 2x (k,u) 



**r J llli J ( 27ri ) 2 
x G (\k + q 2 \,u + rj 2 )G (\k + q x + q 2 \,oj + rj 1 + r] 2 )G (\k + q x \, oo + rji 



x 



(4.1) 



Use of CLHD yields 



T, 2x (k,u) 



it f d 3 qid 3 q 2 



in 



e(|fc + q 1 + q 2 \ - 1)6(1 -\k + 9l |)6(l - \k + q 2 \) 



+ 



0(1 



to — \k 2 + q 1 ■ q 2 — i5 
\k + q t + q 2 \M\k + q 1 \- l)9(|fc + q 2 \ 



u - \k 2 + q x - q 2 + id 



(4.2) 



see also [lOj], Eq. (A. 5). This together with (11.91) used in the Galitskii-Migdal formula (II. lip 
gives (cf. Fig. Id, middle) after the uo— integration has been performed 



v 2x 



3q 4 c 



r Re 



■57T 



d 3 kd 3 q 1 d 3 q 2 



-7 + 



d 3 kd 3 q\d 3 q 2 



Q1Q2 Qi ■ <l2 + W J lid Qi ■ (-Q2) + iS 

B 



• (4-3) 



It is easy to show with the help of the substitutions q 1 — > q[, q 2 — > —q 2 , k — > —{k' + q[) 
that the second term equals the first one. The virial theorem v 2x = 2e 2x gives the 2nd-order 
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exchange energy e 2x to be compared with the 2nd-order direct energy e 2 d [see Eq. f l3.20p ]: 

3g c 4 f d 3 kd 3 qi d 3 q 2 P 3g c 4 f d 3 kd 3 qi d 3 q 2 P 

62X " q \ q \ q x q^ ^ " + V) 3 i qjqj Ql • q 2 ^ 

A A 

P means the Cauchy principle value. Note the replacement l/q\ — > 1/gf and the addition 
of a factor —1/2, when going from the direct term e 2 d to the corresponding exchange term 
e 2x . Having e 2x available, the functions n 2x (k) and Si x ( q ) follow by means of the functional 
derivatives fl!.12j) . 

Combining Eq. ffT~l2l with (jOJ) yields n 2x (k) = j^^F 2x (k) with 

F 2x (k^l) = T ]: [ d3q f 2 q2 ( 1 vj , for |fc + q 1 + q 2 |^l, (4.5) 



(4tt/3) 2 



U) 4 

For comparison the same procedure with e 2 d yields n 2 d(k) = ,. fea- i^d^) with 



F 2d (A; ^ 1) = ±1 / rf3g ^ 2 3g2 1 for |fc + ^ + q 2 | ^ 1, |fc + q lj2 |^l. (4.6) 

^ J <Zi<7i I9i ' Q2) 

It follows from (1C.17I) . (IC.18I) . that (14.61) is equivalent with what arises from (I3.5P for r^2d 
('descreening'). Whereas the direct term drives the electrons outside the Fermi surface and 
decrease the quasi-particle wei ght zp, the exchange process draws them back inside the 



Fermi surface and increases zp 49(. Again the boundary conditions enforce gi j2 — > 00 for 
k — > 00, so the integral becomes simply (4tt / 3) 2 / k 8 . Hence n 2 d(k — > 00) — > +2u; 4 1 /(4fc 8 ) 
and n 2x (k — > 00) — > —cut/ (4k 8 ) rx [n T (k) + n 2x (k)]k-,oo — > +^ 4 1 /(4A; 8 ). Comparison with 
(Oj) shows g(0) « g (0) = 1/2. F 2d , 2x (0) are given in fl022|) . flCT23|) . integral properties of 
F 2x (k) are in ( 1C.24I) . ( 1C.25I) . What concerns N 2x one should expect N 2x « — |JV r , because 
of n 2x (k -> 00) = -|nr(A; -> 00) . Thus JV « iV r /2 + - • • « 0.02691 r 3 / 2 + • • ■ . 

Combining Eqs. (11.121) and (I4.4p . the lst-order exchange term of S(<j) is (Fig. Id, right) 



5ix(g) - + ^j "7"- y^^^ 



(4.7) 



For comparison with S\&{q) the Macke function I(k) in Eq. (13.291) is rewritten with (10. : 



Si d (q) = -2 



^pi y /(g) /(g) _ r d 3 kd 3 q2 P 
4vr/3y q 2 ' q 2 J q \ q l ■ q 2 



(4.8) 
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They have the asymptotics SidJa — > oo) — > —2uj 2 x /q A and S\ x {q — > oo) — > +u 2 x /q A rx 
[S r (q) + <S'i x (q')] g _ >00 = — cjpj/g 4 |40|. Thus the x-term again compensates half of the direct 
term. Comparison with (11. 8ft shows g(0) ~ #o(0) = 1/2, as it should up to this order. How 
does Si x (q) influence the non-analyticity of S(q) at q = 2 and the behavior of S(q) at the 
origin q = 0? 



V. PAIR DENSITY AND SHORT-RANGE CORRELATION 

The SSF S(q) and the PD g(r) are mutually related through the Fourier transforms 

oo oo 

l-*(r) = i d(q 3 ) ^-[1 - S(q)) , 1 - 5(g) = « 3 / d(r 3 ) ^ [1 - g(r)) . (5.1) 
2 J qr J qr 

o o 

5*(0) = expresses the perfect screening sum rule: the normalisation of 1 — g(r) is 97r/4. 
About it, the plasmon sum rule says S(q — > 0) = q 2 /(2uj p i) + • ■ ■ . - For r s = (ideal Fermi 
gas) it is [see also Eq. (12. 3|) ] 

„ , ^ 3 q l/g\ 3 ^/ „ x j /N _ 9 /sin r — r cosr \ 2 _ . . 
5 (g < 2) = -| - - (|J , 5 (g > 2) = 1 ^ (r) = 1 - ^ ( ~ 3 J < 1 -(5.2) 

This causes the potential energy in lowest order as v x = —(3/16) q 2 c = — (3/4) 2 The 
non-analyticity of So(q) at q = 2, namely the 2nd-order-derivative jump ASq(2) = 3/4, 
causes - Fourier transformed - the non-interacting Friedel oscillations of (?o( r )- They are 
tiny: the 1st minimum at r m 5.76 is 1 — 0.0037. 



Short-range correlation means the behavior of the PD q(r) for r l/g c , to which belong 
also (i) the coalescing cusp and curvature theorems 2jJ, |26| and (ii) the influence of the 



on-top PD g(0) on the large- wave-number asymptotics of n(k) and S(q), Eq. ( 11.81) . 



The relations (15. ip between the SSF S(q) and the PD g(r) make, that the on-top value g(0) 
follows from the normalization (1.7) of 1 — S(q). The expansion S(q) = So(q) + S T (q) + 
Si x (q) + ■ ■ • creates corresponding on-top terms gi(0) with i — 0, r, lx, ■ • • . This series starts 
with <7o(0) = 1/2 according to (15.21) . The contribution of S r (q) consists of two parts, a term 
Sid(0), linear in r s , and a term S2v{q), logarithmitically non-analytic. The first order direct 
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term is 



Wo) = 2 / S M 



2 °° 



57T 



-0.73167 r s 



(5.3) 



where ( IC.llj) is used. The factor 2 is killed by gi x (0) = — |<?id(0). So, gi(0) = gid(0) + 
^ix(O) = |(7id(0) ~ —0.3658 r s (32], Exactly this high-density behavior of g(0) results 
also from the ladder theory as a method to treat short-range correlation 5l|, |52J. The second 



term ~ r 2 In r s follows from 



3q f 1 

S 2 r(q) = S T (q) - S ld (q) = — du - ■ 

K J q 



1 R 3 (q,u) 



00 00 



92M = \J d(q 3 ) S 2v (q) = ^ J duj ' qdq 





R 3 (q,u) 



q 2 + q 2 c R(q,u) 



rx 







q 2 + q 2 R(q,u) 



-2 3 



7r 2 \ / 3ar 



4 / V 2vr 



In r s + 



where J du Rl(u) = | (3 — z ^-) is used. Again part (~ hal: 



?) °f 52r(0) is compensated by a 



corresponding exchange term. Thus with x ~ 1/2 it is 32] 

1 , ■ 



9(0) 



a 



7T" 



61n2 - 3)— r s - x 2 ( 3 



7T" 



/3a V 



( — I r s In r s + ■ 



4 y v2tt y 



0.3658 r s - x 0.032966 rj In r s + ■ 



(5.5) 



The decrease of g(0) with increasing r s describes the increase of the area between 1 and 
S(q), the amount of which is 1 — g(0) = 1/2 + 0.3658 r s + ■ ■ ■ . 

VI. SUMMARY 



Following the procedure of Fig. 2, it is shown for the ground state of the high-density 
electron gas (as an example), how the static 2-body quantity S(q), the static structure 
factor (SSF), follows from the dynamic 1-body quantity T,(k,uj), the Dyson self-energy, 
using rigorous theorems as the Migdal formula, the Galitskii-Migdal formula, the generalized 
Hellmann-Feynman theorem, the virial theorem. Along this way all the static RPA results 
are thoroughly revisited and summarized on a unified footing: the energy e [Eq. (13.161) ] and 
its components t and v, the momentum distribution n(k) [Eq. (3.19)], its behavior for k — > 0, 
k — » 00, and at k — > 1 ± with the discontinuity zp, the SSF S(q) [Eq. (3.18)], its behavior 
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at q — > 0, and the on-top pair density g(0). Several identities were found, e.g. the relation 
between the SSF S(q) and the Macke function I(q). S(q) and n(k), stemming from the 
2-body density matrix, are simultaneously linked mutually through the self-energy £(&, u), 
see Fig. 2. An exercise would be to perform the "inverse" procedure S(q) — > v — > e — > n(k). 
So far not solved problems: To have exactly all terms up to ~ r 2 and r 2 In r s available, the 
functions I x (q) and F 2x (k) of Eqs. (1C.19j) and fj4.5j) . respectively, have to be calculated. 
F 2x (/c) has to be renormalized. Besides the exchange term, which compensates part of the 
direct term <?2r(0), Eq. (15.41) . has to be specified and calculated. 
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APPENDIX A: FUNCTIONAL DERIVATIVES 

In terms of diagrams, the relations (11. T 2 [) have a simple interpretation according to 
elementary rules of differentation: Se/Sv(q) means that in a (vacuum) diagram for e one 
after the other interaction line has to be cut out sucsessively: 5v(q')/5v(q) = 5(q — q'). 
So to each e— diagram corresponds a set of S(q)— diagrams. Because t(k) is in the energy 
dominator of Go(k,u) (11.91) . 5e/5t(k) means, that one by one particle-hole line has to be 
cut up in two sucsessively: 5G (k' ,u)/5t(k) = —G (k,uj)S(k — k')Go(k' Thus each 
Go has to be replaced by — GqGq. The minus sign is because a closed loop is broken and 
the 2 Go's are the outer ends of an n{k)— diagram, corresponding to GoSGo, see (13. 4ft . ( 1C.16I) . 

The RPA energy e r of (13.131) is of the kind 

e T [t(k),v(q)] = Jd s qJ dr) f{Q{q,r]),v(q)) (A.l) 



with Q(q,i]) being a functional of t(k). From the derivative 

£e r f df(Q(q r ,r]),v) 
6v(q) 



J dq J d7] *> ^ 5(q - q) = J 



dv 



(A.2) 

v^v(q) 



follows S T (q) of Eq. (I3.14I) . as it should. From the derivative 

6Q{q, V ,[t{k)]) 



~ fd 3 q[d V dmV{q) 



5t(k) J J ' dQ 



*Q(q,v) 



5t(k) 



(A.3) 



follows n T (k) as soon as the task SQ/5t is solved. This is done with the help of 



6 1 --f**, s % + V h) ,r« s ?- h l- (a.4) 



5t(k) J t(k' + q)-t(k')±r]-iS J [t(k' + q) - t(k') ± r] - i6} 2 
This result combined with the definition (11.151) of Q(q, rj) yields 

W-s'**''- (A ' 5) 

This identity used in ( 11. 12ft gives the expression ( 13.151) . which has been derived from the 
RPA-self-energy £ r (fc, a>). 



To derive the identity 

h t(h) 

St{k) 



d3k WttTT I d V Q n+1 (<l,v) = -n I dr] Q n+1 (q,v) ■ (A.6) 
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one must have in mind, that Q(q,r)) has the structure Qf{rj) with f{rj) = l/(r — 77 — iS) + 
1/(t + r] — id) with the excitation energy r = t(k + q) — t(k). After the contour integration 
(only) n factors /j(r n+1 ) remain: 

dr] Q n+1 (q,r]) = Qx ■ ■■Q n+l J dt\ ftirj) ■ ■ ■ f n +i(v) = Qi ■ ■ • Qn+ifi{r n +i) • • -/nfa+i) . 

(A.7) 

Now the operation J c? 3 A; t(k)5/5t(k) counts these remaining n factors. Finally the energy 
denominators l/(Tj + T n +i) cause the minus sign. 

APPENDIX B: THE FUNCTION R(q, u) AND CHARACTERISTIC RPA 

CONSTANTS 

From the particle-hole propagator (11.151) follows the real function R(q,u) = Q(q,iqu) with 

\2 1 „,2 / I , q 1 _ 2' 



1 l+« 2 -^ (| + i)2 + u 2 / 1 + | x _a 

1 H — In -r — u arctan + arctan - 

2q (| — l) 2 + u 2 \ u u 



.(B.l) 



It has the g-expansion R(q, u) = Rq(u) + q R\(u) + ■ ■ ■ with 



1 1 

Ro(u) = 1 - m arctan- , Ri(u) = — — . ... (B.2) 

u 12(1 + u 2 > 2 



Its derivative with respect to u is 

r*, \ 1 / 1 + 1 1 - f \ « , (f + l) 2 + n 2 

ijVqi, u ) = — arctan + arctan H In -§ (B.3) 

2 \ u u J 2q — l) 2 + u 2 

with the g-expansion R'(q, u) = R' (u) + q 2 R[(u) + ■ ■ ■ with 



u 1 . -u . . 32m 

— — --arctan-, ^(u) = — — — ^3, • • • , R {q -> 00, u) = - — 
1 + m z it 3(1 + m 2 ) 13 3g 4 



(B.4) 
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For the integral J du R 2 (q,u) see ( 1C.12j) . Important RPA constants are: 
o 



bl 



7T" 



7T" 



du Rl(u) = - — » 0.031091, 



du Rl(u) \nR (u) » -0.0171202, 



6 r = 6r + a (h---j 



^2d 



3 /" dq 



4?r 4 J g 



-0.045423, 
/(g) 8tt 4 



o 

1 1 

4 7T 2 



9 



-a 6(1 -g) 



11 8 . n 

— --ln2 + 2 In 2 2 
6 3 



-0.025677, 



'2x 



-In 2- « +0.02418. 

6 4 7r 2 



o 

oo 



'/'/ 



1 + W 2 



7T 



—a « 0.482003, / 



Ro(u) ]hRq(u) 

1 +W 2 



^ ^oM arcta n- » -3.353337 
Ro(u) u 





oo 



R"(u) ! 
dw -^-4 arctan - « 4.581817. 
Ro(u) u 



TV" 



1 



1 + U 2 2 



oo 

— I du 



du Ro(u) \nRo(u] 



_4(l + u 



2\2 



'It 



1 

6 



(B.5) 

(B.6) 
(B.7) 



(B.8) 
(B.9) 



Here is a list of integrals: 

oo oo 

/ du i?o( M ) = y a ~ 0.321336, J du R 2 (u)\nR (u) ^ -0.176945, (B.10) 



-0.345751, (B.ll) 



(B.12) 
(B.13) 
(B.14) 



The Kulik function G(a;) approximates F T (k) for A; ~ 1, see (13.421) . It follows from Rq{u) 
according to 

oo 

R' (u) u R {u)-Ro(y) 



G(x) 



du 



Ro(u) u + y u-y 



(B.15) 



y=x 



/y/Ro(u) 



It has the small-x behavior 

G(x < 1) = G(0) + 



7T I ^ + Vs^j x + 0{x 2 ) \nx + 0{x) 



(B.16) 
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with 

oo 

Rn(u) 1 



G(0) = - [ du arctan - « 3.353 337 . (B.17) 

J R {u) u 





The coefficient of xlnx is 7.908 799 (the Kulik number). G(x) has the large x-behavior 

G(x > 1) = -(1 -In 2) — « 0.160 668^r. (B. 

6 x 2 x l 

The Kulik function is shown in (30J, Fig. 1. 

The relative number of particles, N r of (13.471) . follows from Fi(k) of (13.381) as 



oo 



iv r = -M [ du [ qdqk dk R ' M 



4tt J Jj [q 2 + q 2 R(q, u)} 



u u 
arctan — — -- — - — arctan ■ 



o 



(k 2 -l)/(2q) k-q/2_ 

(B.19) 

In lowest order the decisive region is k ^ 1 and < q J$ 5: 



oo 8 1+q 

3 4 f , I , R'(q,u) f , „ 
N, = - T ^]duJ qiq [(j! + qlR(q , a)f J kdk 

1 



u u 
arctan — — ——, — - — arctan ■ 



(k 2 -l)/(2q) k-q/2\ ' 

(B.20) 

The /c-integral yields (expanded in powers of q) 



1 / 1 7T\ 1 , L 1 \ 

- arctan arctan u H — u In 1 H — - o + • • • . 

2 V u 2)2 \ u 2 J 

oo 

The 1st (g-independent) term does not contribute, because of J du R' Q (u)(- ■ •) = 0. The 

o_ 

g-integral yields in the limit r s — > the expression it / (Aq c yJ Rq(u)) , giving ( 13.471) . 

APPENDIX C: IDENTITIES, MACKE FUNCTION, MOMENTUM 

DISTRIBUTION ATfc = 

For the regions of the wave-number integrations the following short hands are introduced 

A = (k < 1, \k + q l + q 2 \ < 1, |fe + <?i )2 | > 1), A = (h, 2 < 1, |fei j2 +q\>l) , (C.l) 
£ = (k > 1, |fc + q l + q 2 | > 1, I* + Qi,2\ < !)> ^ = (A; li2 > 1, |fc li2 + g| < 1) . (C.2) 

For the Heisenberg energy e 2 d, the Onsager energy e 2x , for the SSF's Sid,ix(9) and for the 
momentum distributions n 2 ^ 2x (k) the following equivalent writings exist: 
f d 3 kd 3 q 1 d 3 q 2 P f d 3 qd 3 k 1 d 3 k 2 P 

62d ' J qhl Qi <i2 J q 2 q 2 q- (&i + fc 2 + <?) ' 

A A 
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e2 X : 



d 3 kd 3 qid 3 q 2 P 



d 3 qd 3 kid 3 k 2 



P 



Q1Q2 Q1Q2 J <? 2 0i + k 2 + q) 2 q- (fei + fe 2 + q) ' 

A A 



S ld (q) : / d 3 W 3 g 2 



P 



Qi ■ Q2 



- / d 3 hd 3 k 2 



q- (k 1 + k 2 + q) 



q^q A 



n 2d (k < 1) : 



n 2d (k > 1) : J 



d 3 kd 3 q 2 P 



Q 2 q\ ■ q 2 



d 3 kid 3 k 2 



{k l + k 2 + q) 2 q- (k± + k 2 + q) 



q^q a 



d 3 qid 3 q 2 1 f d 3 qd 3 k 2 P 



(Qi ■ Q2) 2 J <?V [q ■ (ki + k 2 + q)] 2 



d 3 qid 3 q 2 1 f d 3 qd 3 k 



B 



(Qi ■ Q2) 2 J <? 2 <? 2 [q ■ {ki + k 2 + q)] 2 



B 



k\ — >k 



k 1 — > k 



(C.4) 



(C.5) 



(C.6) 



(C.7) 



(C.8) 



n 2x (k < 1) : 



n 2x (k > 1) : 



d 3 qid 3 q 2 1 



d 3 qid 3 q 2 1 



d 3 qd 3 k 2 



P 



q 2 (k 1 + k 2 + q) 2 [q-(k 1 + k 2 + q)] 2 
d 3 qd 3 k 2 P 



k 1 — >k 



(1W2 (9i-g 2 ) 2 J ? 2 (&i + k 2 + q) 2 [q ■ (fei + k 2 + g)p 



B 



(C.9) 



k\ — > k 
(CIO) 



These equivalences are consequences of the invariance, when changing the integration vari- 
ables according to 

k < > ki , q 1 < ► q , q 2 < > -(fei + k 2 + q) . 

Note that Eq. (IC.5I) defines a function I(q), which has been calculated explicitly by Macke: 



7( g ) = - / d 3 kd 3 q 2 



P 



Q1Q2 



in 



or 7(g) = + j d 3 kid 3 k 2 
q^q A 

00 

(4vr/3) 2 



%^0) = — ag+--- , 7(g -> 00) = v ' ' + ••• , / dq 7(g) = — (vr 2 + 6 In 2 - 3) . 



q • (fei + fc 2 + g) 

8tt 2 



45 



(C.ll) 



The linear behavier for g — > causes the Heisenberg divergence of e 2d , the large-g asymp- 
totics makes S(q — > 00) ~ 1/g 4 . At g = 2, 7(g) has an inflexion point and 7"(g) has the 
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jump discontinuity A/"(2) = 2n 2 , besides I(q) = ■ ■ ■ + c±^(q — 2) 3 ln \q — 2| + • • • with 
c + = +1 for q > 2 and c_ = — 1/2 for q < 2 . The non-analyticity at g = 2 influences 
the Friedel oscillations of g(r — > oo). For further properties see [9J, Appendix C, [53J. The 
identity 

oo 

8nq J du R 2 (q,u) = I(q) (C.12) 
o 

is proven as it follows. On the one hand it is (with rj = iqu and R(q, u) = Q(q, iqu)) 



CO 

8tt 2 J ^Q 2 (q,v) = 8Trq J du R 2 {q,u) . 



(C.13) 



On the other hand the lhs of flC.131) - with the definition (11.151) of Q(q,i]) - can be written 
as [with the particle-hole excitation energy n = t(ki + q) — t(ki)] 



(S / M I i 

A 



-T7 + 



1 



T\ — Tj — iS t 2 + rj — \5 T\ + rj — iS r 2 — rj — iS 
The frequency integration is performed as a contour integration with the result 



(C.14) 



8nq I duR 2 (q,u) = I d 3 k 1 d 3 k 2 



[t(kx +q)- t(fcx)] + [t(k 2 + q)- t(k 2 )} ' 



(C.15) 



With t(k) = k 2 /2 the Macke function /(g) of CUB turns out, qed. 



With the help of the functional derivative 

5 1 5(fci + q-fc)-5(fc 1 -fc)+5(fc 2 + q-fc)-5(fc 2 -fc) 



6t(k) q-(k 1 + k 2 + q) 



[q ■ (fei + k 2 + q)] 2 



it is easy to show 
51{q) 



T2 



d 3 ko 



for k h2 ^ 1, |fe 1)2 +q\ ^ 1 



Stikt) J [q-(k 1 + k 2 + q)} 2 

With the help of flC.12p . (IC.13I) . and (IA.5[) it is again easy to show also 

5I(q) 



5t{kt) 



-8qR(q,u)f(k 1 ,q,iqu) . 
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(C.17) 



(C.18) 



( !CH7j) and flCHSl) together with (lC?7jl . flClSl) show the equivalence of fl339l) [n 2d (k) derived 
from the self-energy T> 2d (k,u)} with (14. 6 p [n 2d (A;) derived from e 2d \. 
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Si x (q) : The exchange counterpart of /(g), namely 
d 3 kd 3 q 2 P 



02 Qi ■ Q2 



d 3 k 1 d 3 k 2 



P 



+ k 2 + q) 2 q- (fei + k 2 + q) 



has the asymptotics J x (g — > 00) — > (4-7r/3) 2 /g 2 , the same as I(q). The properties 

00 00 00 

8tt 4 



J dq I x (q) = J dq /(g) and J -f J x (g) 



6 (2tt) 2 
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(C.20) 



follow from 



d 3 kd 3 q 1 d 3 q 2 f d 3 kd 3 qid 3 q 2 



q\ qi -q 2 J q\ Qi ■ <?2 



2 °° 

and \ J dq S lx {q) = v 2x = 2b 2x (ar s ) 2 , (C.21) 



respectively. One may guess, that I x (q — > 0) starts with ~ g 6 . This would contribute to s 4 
ofEq. fl3~33|) . 



n 2d,2x(k) : Here F2d,2x(0) are calculated starting with (j4.5j) . (14. 6p for = 0. The condition 
\Qi < 1 means with ( = cos <(<7x • q 2 ) the restrictions £ = — 1 ■ ■ • — (1 — g 2 — g|)/(2gig 2 ) 
and gi - 1 < g 2 < q% + 1: 



2 gi+l 

_ 2 I / dq 1 



5tt 2 
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w^ijt yWt/* pffSr^— 

,1 1 2 91-I / 

For the corresponding exchange term it is similarly 



(C.22) 
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dqi f dq 2 f dq 1 f dq 2 \ 1 - (g x - g 2 ) 2 



F 2x (0) = +27 r 2 ^ ^+ ^ ^ 



q{ J q$ J q\ 
11 2 



91-1 



02 / 01 + 02 - 1 



-2tt 2 (1 -ln2) 2 



(C.23) 

thus F 2 (0) » -4.1123 + 1.8586 « -2.2537 and n 2 (0) « -0.0063 r 2 . - F 2di2x (A;) have the 
asymptotics — > — 1 1,1 



■ The particle-number conservation 



o / d 3 q 1 d 3 q 2 1 



0i0l (9i ■ 9 2 ) 2 







(C.24) 



easily follows from the replacements in the second integral k — > —k — qi and q 2 — > — q 2 . 
The property 



dk k 



, / d 3 q 1 d 3 q 2 1 
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(C.25) 
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follows from 



oo 



cot r .. ^ e 



(4vr/3) 2 ' d<yk "> F2x ^~ = t2x = _& 2x(ar sy 



(C.26) 



o 



How behaves F 2x (k) near the Fermi edge ? It diverges, but weaker than F 2 d(k), namely log- 
arithmically, F 2x (k — > 1^) ~ ±ln \k — 1|. It is not strong enough to make the corresponding 
(Onsager) energy t 2x divergent: 

- \dkk 4 \n\k-l\ 

1 l-A / 

saying that the lhs of Eq. (IC.26I) exists. The divergence of F 2x (k) is eleminated by the 
renormalization of at least one interaction line. 



~ A In A -> 



(C.27) 



Figures 




FIG. la: Feynman diagrams for the one-body Green's function of the ideal Fermi gas 
Go(k,u), the bare Coulomb repulsion Vo(q) = 7f 2 q 2 /q 2 , and the RPA polarization propa- 
gator Q(q,rj) as defined in Eqs. ( I1.13I) - (I1.15I) . 




FIG. lb: Feynman diagrams for v r — Vq + voQv r , the screened Coulomb repulsion in RPA. 






FIG.lc: Feynman diagrams for the ring-diagram-summed self-energy S r (left) with n r = 
GoS r Go, for the energies e r , v T (middle), and for S r (right) according to Eqs. ( 13. II) . ( 13.41) . 
( I3.13p . (I3.14p . "Descreening" v r — > vq makes the divergent quantities S 2 d, n 2 d, ^2d, v 2 d, Sid- 
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FIG. Id: Feynman diagrams for S 2x (left) with n 2x = Go^xGo, for the energies e2 X , ^2x 
(middle), and for Si x (right) according to Eqs. (I4.1l) - (l4.5p . (14. 7L 




FIG. 2: Terms following from the off-shell self-energy £(/c, lo) as a functional of t(k) = k 2 /2 
and v(q) = q 2 /q 2 '- The Migdal formula (11.101) yields the momentum distribution n(k), from 
which follows the kinetic energy t. The Galitskii-Migdal formula (II. lip yields the potential 
energy v. The Hellmann- Feynman functional derivative ( 11. 12ft of the total energy e = t + v 
yields the SSF S(q), its Fourier transform gives the PD g{r). The cases E x , S r , S 2x are 
studied in Sees. II, III, IV, respectively. 
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FIG. 3: The momentum distribution (13.381) and its /c-g-plane. 



